Mechanisms of degeneracy breaking in pyrochlore antiferromagnets 
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Motivated by the low temperature magnetization curves of several spinel chromites, we theoret- 
ically study classical mechanisms of degeneracy lifting in pyrochlore antiferromagnets. Our main 
focus is on the coupling of spin exchange to lattice distortions. Prior work by Penc et al. (Phys. 
Rev. Lett. 93, 197203 (2004)) has demonstrated that such coupling leads to a robust magnetization 
plateau at half the saturation moment per spin, in agreement with experiment. We show that a 
simple Einstein model incorporating local site distortions generates a "universal" magnetic order 
on the plateau, and highlight the distinct predictions of this model from that in Penc et al. (Phys. 
Rev. Lett. 93, 197203 (2004)). We also consider the complementary degeneracy-lifting effects of 
further neighbor exchange interactions. We discuss the implications for transitions off the plateau 
at both the high field and low field end, as well as at fields close to the saturation value. We predict 
that under certain circumstances there is spontaneous uniform XY magnetization (transverse to the 
field) for field values just above the plateau. These features may be tested in experiments. While 
selecting a unique magnetic order in the half magnetization plateau, at zero magnetic field the Ein- 
stein model retains an extensive degeneracy, though significantly reduced compared with the pure 
Heisenberg antiferromagnet. 

PACS numbers: 75.10.-b,75.10.Jm,75.25.+z 



I. INTRODUCTION 

The pyrochlore lattice with nearest-neighbor anti- 
ferromagnetically coupled spins is well-known as one of 
the most frustrated and degenerate magnetic systems^*^ 
Ultimately, this degeneracy must be lifted at low temper- 
ature, but the mechanisms responsible can vary greatly 
from material to material and also depend on applied 
fields, pressure, and other variables. In this paper we fo- 
cus primarily on how degeneracies present at finite mag- 
netic fields are lifted by the coupling of spin and lattice 
degrees of freedom. 

Recent experiments on a number of insulating chromitc 
compounds, namely ZnCr204, CdC^CU and HgCr204 
have revealed distinctive common features in their low- 
temperature magnetization curves^ and other interest- 
ing properties in neutron scattering!^ At low magnetic 
fields the magnetization curve grows linearly with mag- 
netic field. At one point there is a sharp jump in mag- 
netization onto a rather robust plateau, with half the 
full magnetization per spin. In HgCr 2 04 it is possible 
at yet higher fields to observe a smooth transition off of 
the half-magnetization plateau, and a gradual increase 
in magnetization up to what may be a fully polarized 
plateau stated 

The Cr +3 ions sit at the center of octahedra of 0~ 2 
ions, and thus the outer d-orbital electron shell under- 
goes crystal field splitting to a lower energy tig orbital 
triplet, and a higher energy e g orbital doublet. The ti g 
orbitals hold 3 electrons, and therefore by Hund's rule 
form a spin | degree of freedom, with no orbital degener- 
acy (therefore the cooperative Jahn-Teller effect cannot 
lift the degeneracy in this system). These spins are the 
source of magnetic behavior in these compounds. The 



Cr +3 ions sit on the sites of a pyrochlore lattice, and 
therefore a minimal model for the magnetic properties 
of these compounds is the nearest-neighbor Heisenberg 
anti-ferromagnet, with the Hamiltonian 



n 



H 



(i) 



Here H is proportional to the applied magnetic (Zeeman) 
field. In this paper, we will treat the spins as classical, 
and simplify notation by normalizing them as unit vec- 
tors (absorbing a factor of S 2 into J). Much previous 
work has been devoted to this and similar models in zero 
magnetic field &2*i2iiiii2ii&±i A useful rewriting of the 
Hamiltonian is 



W=^]T[(St-h) 2 -h 2 ] 



(2) 



where St = Yljet &j 1S the sum of spins on a tetrahedron 
labeled by t, h = H/2J = hz and we have ignored a 
trivial constant term in the Hamiltonian. This model 
has a macroscopically degenerate classical ground state 
manifold at all fields up to full polarization: any state 
with St = h on all tetrahedra is a classical ground state. 

Within this (over-)simplified picture, the magnetiza- 
tion is everywhere a smooth (linear at low temperatures, 
since the magnetization is proportional to the average St) 
function of the field, and there are no plateaus. Instead, 
some other effects or interactions must be considered to 
explain the observed plateau. On general grounds, a 
plateau is expected to be associated with some break- 
ing of degeneracy, into a state in which the spontaneous 
static moments of the spins are all parallel or antiparallel 
to the field axis. Quantum mechanically, this is simply 
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FIG. 1: The pyrochlore lattice. A network of corner sharing 
tetrahedra. 



because, in a non-collinear ordered state, the Goldstone 
theorem associated with U(l) symmetry-breaking around 
the applied field direction ensures the presence of a gap- 
less magnetic excitation (spin wave), which contradicts 
the incompressibility of the magnetization plateau. The 
analogous classical argument is that a non-collinear spin 
state may always be arbitrarily slightly deformed to lower 
its free energy in response to a change in field, which 
by thermodynamics implies a non-constant magnetiza- 
tion. Hence, to understand the plateau, we seek mech- 
anisms to select one or a set of collinear ground states 
out of the classically degenerate manifold. On the py- 
rochlore lattice, the natural collinear states for the half- 
polarized plateau£*i£*i& are those with 3 "majority" spins 
aligned parallel to the magnetic field, and 1 "minority" 
spin aligned antiparallel to the field, on each tetrahe- 
dron. However, even if one assumes a collinear state for 
the spins, a massive degeneracy of the ground state still 
remains, since there is considerable freedom in fixing the 
location of the minority spin on each tetrahedron. 

The possibility that quantum fluctuations might con- 
trol the state selection - of and within the collinear 3:1 
manifold - has been explored elsewhere i^i^iL Here we 
will investigate alternative mechanisms, within classical 
models^ A guide to the possible physical processes in- 
volved comes from two sets of observations. First, it has 
been noted experimentally that the above chromite ma- 
terials, ACr204, exhibit strong magnetostriction, espe- 
cially upon entering the plateau region4£ This strongly 
suggests that spin-lattice coupling plays an important 
role in the plateau formation. Second, studies of the 
structurally and electronically analogous set of spinels, 
ACr2(S,Se)4 - with S or Se replacing O atoms and the 
same non-magnetic A atoms - display ferromagnetic ten- 
dencies or long-range order, and in some cases an appar- 
ent competition of ferromagnetic and antiferromagnetic 
interactions^ This indicates drastic changes in the mag- 
netic interactions may be effected by small changes in 



structural parameters. More specific implications of the 
trends in these materials for the chromites will be dis- 
cussed later. 

This main analysis and results of this paper are as 
follows. Guided by the above observations, we focus pri- 
marily on a minimal model for the plateau structure in- 
volving only spin-lattice coupling. In this minimal model, 
the lattice modes are taken into account by the simplest 
possible Einstein phonons describing motions of the mag- 
netic sites. We show that this model indeed captures a 
simple and robust mechanism for plateau formation and 
predicts a unique ordered 3:1 state - the "R" state, shown 
in Fig01 - on the plateau. Extended to the full range 
of magnetic fields, this Einstein model predicts a first- 
order transition to a non-collinear ground state at lower 
fields below the plateau, and a second-order transition 
to a canted ferrimagnetic state above the plateau. The 
canted ferrimagnet retains the Ising order of the R state, 
but in addition possesses XY ferrimagnetic order of the 
magnetic moments transverse to the field axis. At zero 
field, the Einstein model retains a large ground state de- 
generacy, though it is still vastly reduced from that of 
the ideal model without spin-lattice coupling. 

A plateau with the same R state structure can also be 
stabilized by a combination of spin-lattice and further- 
neighbor exchange interactions. We give the conditions 
on these further-neighbor exchanges for this to occur. 
Consistency of this more complicated but still feasible 
scenario could be then tested by placing independent 
constraints on these couplings from other measurements. 
This is considered further in the Discussion. 

A number of studies of spin-lattice and further- 
neighbor exchange effects in pyrochlores have already 
been carried out. A well-known analysis of certain zero 
field spin-lattice couplings by Tchernyshyov et al. chris- 
tened the resulting degeneracy-breaking a "spin Jahn 
Teller" effect iSS Because this analysis was at zero field, 
and because it considered only q = phonons, it has 
little bearing on the present work. More relevant is the 
pioneering study of spin-lattice couplings on the plateau 
by Penc et aii Their work provides a simple explanation 
of the plateau formation, but unlike the theory in this pa- 
per, does not explain the breaking of degeneracy within 
the 3:1 states. We will compare their "bond phonon" 
Hamiltonian with our Einstein model throughout this pa- 
per. 

The remainder of this paper is organized as follows. In 
Sec. [n] we present two models of spin-lattice coupling, a 
"bond" model and a "site" model. In Sec. IHII we discuss 
the implications of these two models for the magnetic or- 
der on the half-saturation magnetization plateau and the 
transition off the high and low field edges of the plateau. 
We discuss a more general model with further neighbor 
spin interactions in Sec. IIVI Finally, a discussion of our 
main results and their relevance to experiment is given 
in Sec. E| 



3 



II. SPIN-LATTICE COUPLING 

In this section we discuss some simple models for the 
coupling of the magnetic degrees of freedom to phonon 
modes. We will treat the spins and phonons classically 
and in equilibrium. With these assumptions, the statis- 
tical mechanics of the phonons is captured by a Gaus- 
sian integral over the associated set of displacement co- 
ordinates in the partition function. In such a case, the 
phonons can (if desired) be integrated out to obtain an 
effective spin Hamiltonian which contains additional in- 
teractions beyond the Heisenberg form. 

Let us first make a few general comments regarding 
spin-lattice interactions. For a fixed, static, distortion of 
the lattice, we expect modifications of the exchange in- 
teractions that are (to an excellent approximation) linear 
in the displacement coordinates. Neglecting weak spin- 
orbit effects, the exchanges remain of Heisenberg form. 18 
Therefore the general form of the modified exchange is 



(ij) 



Sj [1 - -yuij] , 



(3) 



where «y is the linear combination of displacement co- 
ordinates coupled to the pair of spins i,j. Because all 
nearest-neighbor pairs of pyrochlore sites are equivalent, 
they are all described by a single spin-lattice constant 7. 

A naive interpretation of Eq. © is that tiy is pro- 
portional to the distance between spins i and j. As the 
distance is increased/decreased, the overlap between elec- 
tronic wavefunctions on the two sites decreases/increases, 
leading to a change in the exchange coupling, propor- 
tional to this distance. This picture is in fact appropri- 
ate for direct exchange, in which there is no intervening 
oxygen as in superexchange. In the spinel chromites, the 
antiferromagnetic Cr-Cr is indeed believed to arise from 
direct exchanged More generally, the dependence of ex- 
change on displacements may involve changes in the bond 
angles as well as distances. Nevertheless, from this sim- 
plistic view, one expects 7 > (hence the minus sign in 
Eq.©), and 7 of order the inverse of the effective Bohr 
radius of the electronic orbitals involved. 



A. Bond phonon model 

In the model of Penc et al^^ the Uij £1X6 taken as in- 
dependent parameters, i.e. the length of each pyrochlore 
bond can independently expand or contract. This "bond 
phonon" (BP) model has the elastic energy 



T-/BP _ 



ks_p \ - 

2 

(ij) 



(4) 



where fcgp is an elastic constant. Because each cou- 
ples only to a single nearest-neighbor pair of spins, and in 
this model each Uij is independent, the spin-lattice inter- 
action does not induce any effective interactions amongst 



further neighbor spins. Instead, integrating out the «,j 
according to 



e -0«eff [{Si}] 



(ij) 



. (5) 

one obtains, up to a constant, an effective spin Hamilto- 
nian of the form 



BP 
eff 



./^'S, -S, /;:S, -S.r 



(ij) 



(6) 



Thus, in this BP model, the spin-lattice coupling induces 
an effective bi-quadratic interaction or relative strength 

2 T 

b = ^ — between nearest-neighbor spins. Note that this 
term favors configurations in which neighboring spins are 
either parallel or antiparallel, i.e. collinear configura- 
tions. The BP model therefore gives a simple explanation 
for the preference for 3:1 states in the field range in which 
the classical Heisenberg model prefers half-magnetization 
states -liiii 

The preference for collinear spin arrangements can be 
understood physically in terms of the phonons as follows. 
If a given pair of spins is antiferromagnetically aligned, 
then the spin-lattice coupling in Eq. © can be made most 
negative by choosing < 0, i.e. contracting the bond 
to enhance the effective exchange. Conversely, if a pair 
of spins is ferromagnetically aligned, the bond can ex- 
pand (u^ > 0) to weaken the ferromagnetic exchange 
interaction. In either case, the bond energy is lowered by 
the same amount (because of the linear phonon coupling) 
relative to the undistorted bond. 

It is straightforward to see that, as claimed earlier, all 
the 3:1 plateau configurations remain degenerate within 
the BP model. By rewriting the exchange interaction as 
in Eq. (|2J), one obtains 



J 



L t 



h) 2 -h 2 ]-6J^(S i .S,) 2 . (7) 

(ij) 



For h — 2, every 3:1 configuration minimizes the first 
term on each tetrahedron as well as the second term on 
each link. Hence they are the global ground states and 
all degenerate. Furthermore, even for different values of 
h, all these states remain degenerate, since they have the 
same (no longer minimal in exchange energy) St on each 
tetrahedron, and the same (minimum in energy) value 
of the bi-quadratic term. We will return to discuss the 
magnetization curve and configurations away from the 
plateau in Sec IIIII 



B. Einstein (site) phonon model 

The lack of splitting of the degeneracy of the 3:1 states 
is a non-generic feature of the BP model. It arises from 
the fact that the bond displacements are taken to be 
completely independent of one another, so that they can 
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induce no spin correlations beyond nearest-neighbor. In 
reality, however, this is not the case. To make a change 
in a given bond length requires moving one or both of 
the atoms involved, which will at least distort the other 
bonds connected to these atoms. A more natural phonon 
model can be formulated in terms of the independent dis- 
placements of each atom, with the bond distances deter- 
mined from these atomic displacements. If the harmonic 
phonon energy is taken to be a sum of independent restor- 
ing forces for each atom, this is simply the conventional 
Einstein model. As usual, such an Einstein model pro- 
vides a crude but reasonable approximation, provided the 
most important phonons are optical phonons rather than 
the long-wavelength q w acoustic modes. 

We therefore adopt this Einstein (site rather than 
bond) phonon model. To derive an appropriate form, let 
us assume a distance dependent exchange coupling (see 
e.g. Ref. 1 1 SI in a zero field context). It can be expanded 
in small atomic displacements Uj for each site i: 

Jit = J(\ri - Tj \) « J(R) + (m - Uf) ■ VJ(R) + . . . 
«(l-7eij • (iij — Uj)) J , 

(8) 

where R is the vector between the pyrochlore sites i and 
j, and &ij — R/|R| is the unit vector in the corresponding 
equilibrium direction. Comparing to Eq.Q, the Einstein 
model has then 

Uij = Gij ■ (Ui - Uj) , (9) 

with the elastic energy given by 

i 

It is amusing to note that although the independent 
modes considered are quite distinct in the bond and site 
phonon models, on the pyrochlore lattice both formu- 
lations contain the same number of degrees of freedom. 
The site displacement vectors represent 3 degrees of free- 
dom for each site, which for N sites of the pyrochlore 
lattice gives a total of 3N degrees of freedom. On the 
other hand, there are 6 links per tetrahedra, and one 
scalar bond phonon per link. There are N/2 tetrahedra 
in the lattice, so the total number of degrees of freedom 
in the bond phonon model is also 3N. 

Just as we did for the BP model, we now proceed to 
integrate out the site Einstein phonons. Because the 
Hamiltonian is quadratic in the displacements, this Gaus- 
sian integration is equivalent to simply minimizing the 
Hamiltonian with respect to the set of Ui and eliminating 
these in favor of the spin variables. The optimal values 
of the displacements are simply 

^-T 2 £ ( S *-S;)e y -> (11) 

kE ieN(j) 

where N(j) denotes the set of nearest neighbors of site j. 
At this point, it is evident already that the presence of 



lattice distortions is tied to frustration. If all exchanges 
could be satisfied equally, i.e. (Sj • Sj) = const, then the 
distortion vanishes: (u*) ~ J2ieN(j) e *i = 0. A distorted 
lattice is thus induced only in frustrated states, and for 
instance no distortion is expected at large fields where 
the spins are fully and uniformly polarized. If a distortion 
is still observed in a fully polarized state, it cannot arise 
from this mechanism, and might have no connection with 
the magnetic order of the system. 

Substituting back Ea. (|ll|) in the Hamiltonian, we ob- 
tain the effective Hamiltonian for the (site) Einstein (E) 
model: 

Wf ff = j£S s ;. Sj -^]r|u*| 2 , (12) 

(ij) 3 

where, as above, the Zeeman interaction with the ex- 
ternal field has been dropped for brevity. The form in 
Eq. (|12|) is actually quite convenient for analysis, but we 
first write out the induced interactions explicitly for com- 
parison with the BP model. One finds 

Hf ff = JJ2 [Si-Sj - b (S.-S,) 2 " 

(ij) 
b' 

~ J ~n ^ ] (Si ■ Sj)(S; ■ Sfc)ejj • Bjfej 

(13) 

with b = as before. Here we have defined a sepa- 
rate parameter b' which, according to the strict devel- 
opment above, is not independent (£>' = b). However, 
it distinguishes the additional terms in the effective spin 
Hamiltonian which are not present in the BP model. We 
propose to view b' as a separate phenomenological param- 
eter to describe the phonon modes. In particular, taking 
< b' < b corresponds to "softening" the bond phonons 
somewhat, and interpolating between the Einstein model 
and the BP model. 

Because the same biquadratic term is present as in the 
BP model, we may expect that, as in that case, collinear 
states are favored. Indeed, when the additional terms are 
weak, b' <C b, they may be viewed as weak perturbations. 
Since the b term already selects collinear 3:1 configura- 
tions, the b 1 > term can then only select a ground 
state within this set of states, and has little effect on 
the stability of the plateau. We show below that b' > 
selects a particular collinear classical spin configuration 
depicted in Fig. 4, in the manifold of 3:1 states. Our anal- 
ysis shows that for the range of parameters b'(< |) the 
plateau retains a finite width, and the same particular 
3:1 configuration prevails. 

To get some intuition for how additional degeneracy- 
breaking is induced by the b' terms, assume a collinear 
state, Si = zoi, with <Ji = ±1 Ising spins satisfying the 
3:1 constraint. The first line in Ea. becomes con- 
stant, and the four-spin product (Si • S,-)(Sj • Sfc) = <JjO~k 
becomes an effective two-spin interaction. Then the 
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(a) Site adjoining the 
tetrahedra is the only 
minority site. 



(b) Two minority site 
connected by two 
parallel links. 




(c) Two minority site 
connected by two links 
bending. 



(d) Distortion of the 
bending configuration. 



FIG. 2: (Color online) The three generic configurations of 
minority sites (red arrow) on two adjacent tetrahedra. 



Hamiltonian within the 3:1 manifold takes the form of 
Ising exchange terms: 



H 



3:1 



X 



off 



am 



X 



cir 



«<«>» 



(14) 



where we have dropped some terms which are constant 
in the 3:1 manifold. The parameters J$ s = Jb'/4 and 
J| ff = Jb'/2 are the effective Ising exchanges between 
second-neighbor sites (connected by two bent links as in 
Fig |2(c)| ) and third -neigh bor sites (connected by two par- 
allel links as in Fig 2(b) and explicitly shown in Fig. [31). 
Note that J| ff is twice as large as J| ff , so the third neigh- 
bor interactions tend to be favored over second neighbor 
interactions in the determination of the magnetic state. 

Clearly these effective interactions break the degener- 
acy of the 3:1 manifold. To understand this breaking 
more physically and thereby derive the selected ordered 
state on the plateau, we return to the formulation of 
Eci. (|12[) . From Eq. l|12l) it is readily apparent that the 
larger the induced lattice displacement, the lower the en- 
ergy. For 6 < 1, we may assume half-magnetization on 
each tetrahedron in the plateau region, and then consider 
the effect of the lattice displacement as a perturbation. 

Let us begin by first assuming the 3:1 constraint on 
each tetrahedron. We then wish to understand which 




FIG. 3: (Color Online) Further neighbor interactions in the 
pyrochlore lattice. The nearest neighbor interaction Ji is be- 
tween site a and site b, the next nearest neighbor interaction 
J-2 is between site a and site c, and finally the next next near- 
est neighbor interaction J3 is between site a and site d. 



configurations of minority spins maximize the displace- 
ments. Consider 2 tetrahedra adjoined by a site j. There 
are 3 generic configurations of the positions of the minor- 
ity sites, depicted in Fig. [3 A fixed fract ion (1 /4) of all 

therefore 



configurations are necessarily those in Fig 2(a) 



the energy of this configuration is irrelevant to the split- 
ting. Of the remaining two configurations, it is simple to 
understand that a nonzero u : 
figuration Fig 



2(c) 



. 3 can only arise in the con- 
because it is the only configuration 
with asymmetry about the site j. Therefore, we wish to 
maximize the number of configurations of this type. 

We can describe this favored configuration with a 
"bending rule". For every pair of adjacent tetrahedra 
adjoined by a majority site, mark the links connecting 
between the two minority sites. These marked links form 
paths on the lattice, connecting all the minority sites. It 
is energetically preferable for these these paths to bend, 
rather than continue on a straight line. Clearly, the max- 
imum number of such bent paths is obtained if all paths 
are bent, i.e. all minority sites are in the configurations 
in Fig 2(c) This can indeed be achieved. By careful 
enumeration of all configurations ( see appendix [5] f° r 
details), it can be shown that there is a unique (up to lat- 
tice symmetries) configuration which satisfies this "bend- 
ing rule" on a "pyrochlore cell" - a volume containing 4 
hexagonal plaquettes (see Fig.^J. If this configuration is 
required on all such cells, the ground state is completely 
specified. This is precisely the "R" state obtained in 
Ref. in a very different quantum dimer model analy- 
sis. The "R" state has space group P4332 and may be 
thought of in terms of filling the pyrochlore lattice with 
a fraction 1/4 of hexagonal plaquettes with alternating 
up/down spins and a fraction 3/4 of plaquettes with one 
down spin and five up spins. 

Because the argument above proves the R-state is the 
best possible collinear state, any alternative ground state 
must be non-collinear. Since non-collincar states cannot 
exhibit a plateau, proving that the R-state is the global 
minimum energy state is equivalent to proving the exis- 
tence of a plateau. However, it is important to emphasize 
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FIG. 4: (Color online) Spin configuration of the R state. Ma- 
jority sites are colored yellow, minority sites are colored blue. 
The flippable plaquette in this unit cell of the R-state config- 
uration is highlighted in red. 



that the above argument assumes the 3:1 configurations 
(actually it can be made equally rigorous assuming only 
collinearity). For b' <C b, this assumption is valid, and the 
above argument becomes controlled. For the generic sit- 
uation with b' ~ &, the effective Hamiltonian is actually 
"frustrated" in the following sense. Because the Einstein 
phonon displacement resides on a pyrochlorc site and is 
related via Eq. Qllfl to spins on the two neighboring tetra- 
hedra, the natural unit for the effective Hamiltonian is 
no longer a bond but such a pair of adjacent tetrahedra. 
One may rewrite the Hamiltonian as a sum over such 
pairs, parametrized by the pyrochlore site j they share: 



ith 



J r 



(s tl - h) 2 - h 2 



J r 



(St, - h) 2 - h 2 



(ife)e*i,2 



(15) 



J- ( S » ' S j)( S j • s k)e jl • e jk 



The extra factors of j above relative to Eq.(J2Jl in the 
nearest neighbor terms are due to the fact that every 
tetrahedron is shared by 4 pyrochlore sites. One can 
show that for b, b' > 0, for a non-zero window of fields 
h in the neighborhood of h = 2 (the plateau region), Hj 



is minimized, when the pair of adjacent tetrahedra shar- 
ing site j is in the "bending" configuration pictured in 
Fig 2(c) Clearly, however, this condition cannot be si- 



multaneously satisfied on every tetrahedral pair, because 
some of the tetrahedron pairs must have a minority site 
adjoining them. Thus the Einstein site phonon model 
exhibits "tetrahedral-pair frustration". The R-state ar- 
gued for above resolves this frustration in a natural way, 
by minimizing the energy on a maximal fraction of tetra- 
hedral pairs (which is 3/4). Because such a relatively 
large fraction of tetrahedron pair units are "satisfied" , it 
appears plausible that the R-state is indeed the global 
ground state. 

We have searched numerically to check for the alter- 
native possibility, that a lesser fraction (perhaps zero) of 
units are satisfied, but that the energy of the dis-satisfied 
units is sufficiently better as to lower the total energy. In 
order to take the "tetrahedral-pair frustration" into ac- 
count, we must go beyond the above single-tetrahedron 
analysis for the BP model. In particular, we must con- 
sider units larger than a single tetrahedron, and also 
larger than a single tetrahedron-pair unit: since these 
units are frustrated, they cannot be simultaneously sat- 
isfied at most fields. Instead, we consider a cluster of 
five tetrahedra consisting of a central tetrahedron and 
its four surrounding neighbors. This is the smallest col- 
lection of tetrahedra for the 3:1 stat es in which a pair 
in the unsatisfied configuration of Fig. 2(a) must appear. 
The Hamiltonian can be written as a sum over the up 
pointing tetrahedra TL — Etet^*' * being the central 
tetrahedron in each cluster. The down pointing tetrahe- 
dra are counted in 4 different clusters in this scheme, so 
we account for this by writing the cluster Hamiltonian as 



Hi 



J r 

'2 



(Sto-h) 2 -^' 
(S t2 - h) 2 - h 2 



J 



(S tl - h) 2 - h 2 
(St, ~ h) 2 - h 2 



V 



(S t4 - h) 2 - h 2 \ -J- J2 ( s * • 

(ife>Sto,l,2,3,4 

4 



J ~9 Y Y ( Sl ' S J')( S j • S fc) e ji • e 3k , 



(16) 



where to is the middle tetrahedron, tx..4 are the other 
four tetrahedra, and the sites j — 1...4 are the 4 corners 
of the tetrahedron to. 

The conclusion from our numerical study, is that for 
< b' < 6/2, the minimum of Ti t above is indeed a 
3:1 configuration comprising a corresponding section of 
the R-state. Hence, because such a configuration can be 
simultaneously realized on every such unit, in this pa- 
rameter range, the spin-lattice coupling indeed stabilizes 
a state with the R-state symmetry. The width of the cor- 
responding plateau is discussed in the following section. 
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III. AWAY FROM HALF-POLARIZATION 

In this section, we explore the properties away from 
the magnetization plateau. 

A. BP Model 

Let us first review the findings of Penc et aim in the 
BP model. The basic results can be understood by simple 
considerations on a single tetrahedron. Such a simplifica- 
tion is satisfactory in this case because the Hamiltonian 
can be written as a sum of such terms on each tetrahe- 
dron, and they can be simultaneously satisfied. Thus the 
BP model does not suffer from "tetrahedral-pair frustra- 
tion" . 

For magnetization greater than half polarization, the 
ground state has a 3:1 configuration, with 3 majority 
spins and 1 minority spin on each tetrahedron. How- 
ever, they are not collinear, except on the plateau and at 
saturation. These vary in such a way that the net spin 
per tetrahedron is increased from 2, by smoothly rotating 
the minority spin from down to up. Because the collinear 
plateau state can be smoothly deformed into the states 
above the plateau, there is a continuous transition at the 
upper plateau edgei All these states above the plateau 
share the same degeneracy as the plateau states: the lo- 
cation of the minority spins is not determined in the BP 
model. 

On the other hand, a state with magnetization per 
tetrahedron of less than 2 cannot be achieved with a 3:1 
configuration. Thus, for fields below the plateau, the 
structure of the configurations must change. Instead, 
over most of this region of phase space, the spin configu- 
ration assumes a 2:2 form. This implies a discontinuous 
change of spins, and gives rise to a first order transition. 1 
Like the 3:1 states, these 2:2 states are highly degenerate, 
due to the many equivalent manners in which each of the 
2 spin polarizations may be arranged. Both above and 
below the plateau, the low temperature phases break ro- 
tational symmetry about the field axis, but with no net 
moment in the x — y plane. 
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FIG. 5: (Color online) Magnetization curves for various values 
of b' with 6 = 0.1 fixed, where b and b' are given below Eq. 1151 . 
The plateau width decreases with growing b' . 



the Einstein model, however, is much more involved than 
the above single-tetrahedron analysis for the BP model, 
due to the "tetrahedral-pair frustration" , explained in the 
previous section. Following the above treatment, we con- 
sider a collection of five tetrahedra. On this cluster, we 
numerically minimize the energy (|l(jf) for each field, and 
determine the zero-temperature M(H) curve for given 
values of b, b' . The magnetization curves for one value of 
b = 0.1 and various values of b' are plotted in Fig. EI 

In the magnetization curve shown in Fig. [3] an abrupt 
jump is observed going from the linear regime into the 
half magnetization plateau region, in accordance with 
the results of Penc et ali Given the 2:2 spin configu- 
ration at zero magnetic field, and the 3:1 configuration 
in the plateau region, we can contemplate a sharp transi- 
tion onto the plateau between 2:2 and 3:1 configurations. 
Such a transition would have to be first order, as the 
symmetry groups are not Landau-Ginzburg compatible 
for a second order transition (one is not the subgroup of 
the other). 



2. Phases 



B. Einstein model 

1. Magnetization 

The BP model captures rather well the broad behavior 
of the low-temperature magnetization curve, M(H), in 
the chromites, where it has been observed. The only 
qualitative exception is the observation of a small feature 
at H m 37T for HgCr204, in the field range between the 
plateau and saturation, which has been suggested as an 
additional phase transition^ 

The magnetization curves we obtained for the Einstein 
model (b' < b) are similar to those for the BP model over 
most field values. Studying the magnetization curve for 



As we have seen above, the differences in the magne- 
tization curves of the Einstein and BP models are min- 
imal. The magnetization, as a purely thermodynamic 
quantity, is only weakly sensitive to the detailed correla- 
tions between spins beyond one or two lattice spacings. 
A better test of the differences between the models is 
to compare their phase diagrams. In the BP model, as 
described above, their are four "phases" , in which the 
local structures on each tetrahedra are distinct: at low 
fields, a 2:2 structure or a "splayed" structure; on the 
plateau, a collinear 3:1 structure; above the plateau, a 
non-collinear 3:1 structure, and at high fields, the fer- 
romagnetic fully saturated configuration^ We have used 
quotation marks around the word "phases" because all 
but the ferromagnetic configuration exhibit an unphysi- 
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Number of unit cells 


Number of bending ice configurations 


2x2x2=8 


12 


2 x 2 x 4 = 16 


36 


2 x 4 x 4 = 32 


82 


4 x 4 x 4 = 64 


216 



TABLE I: Bending Ice configurations. Dimensions of py- 
rochlore cluster indicated. Periodic boundary conditions were 
used. 



Number of unit cells 


Number of ice configurations 


2x2x1=4 


78 


3x2x1=6 


534 


2x2x2=8 


2970 


3x3x1=9 


7974 


5 x 2 x 1 = 10 


28326 


3 x 2 x 2 = 12 


87684 



TABLE II: Ice configurations. Dimensions of pyrochlore clus- 
ter indicated. Periodic boundary conditions were used. 



cal macroscopic degeneracy not related to symmetry. 

A full determination of the phase diagram in the Ein- 
stein model is beyond the scope of this paper. However, 
we will outline those features which are similar and those 
which are clearly distinguishable from the BP model. At 
zero field, it can be shown that the ground states of the 
Einstein model are far less degenerate for all b' > than 
those of the BP model. To see this, we use the repre- 
sentation in Eq. l|15fl . and consider the minimum energy 
configuration on a single tetrahedral-pair with h = 0. 
Simple analysis shows that this minimum energy occurs 
for collinear states with a 2:2 ratio of "up" and "down" 
spins on each of the two tetrahedra, satisfying the "bend- 
ing rule" if links are drawn between the spins aligned with 
the central one. Because the field h = 0, the spin axis 
is arbitrary. Unlike in the plateau region, such 2:2 states 
are unfrustrated: every tetrahedral pair can be chosen 
to have such a configuration. In fact, these states are 
still highly degenerate. They correspond to "ice-rules" 
states, with the additional constraint of the "bending 
rule" . Though we do not have an analytical calculation of 
the number of such states, we have performed a numer- 
ical enumeration of them on various finite size clusters 
(see Table [fl . Evidently, the number of such "bending 
ice" states grows rapidly with system size. It is likely 
that these states are macroscopically degenerate. Nev- 
ertheless, this set of states is much less degenerate than 
the zero field ground states in the BP model, which are 
simply all the 2:2 "ice" states, without the bending rule 
imposed (see Table HTjl . We will comment upon the phys- 
ical consequences of this degeneracy in Sec. [V] 

At small non-zero fields, we expect the same "bending 
rule" states to remain approximate ground states, with 
the axis of the two spin orientations "flopped" into a fixed 
one at a small angle (proportional to h/ J) away from the 



x — y plane in spin space. Thus in this region there is a 
broken XY symmetry around the spin-rotation axis. In- 
deed numerical minimization of a tetrahedron pair shows 
that for finite weak magnetic field the nearly collinear 
2:2 bending state persists. For intermediate fields half- 
way between zero field and the plateau, we do not have 
definitive results. As in the BP model, there is a first 
order transition separating the low-field region from the 
plateau. 

Like in the BP model, at fields just above the plateau, 
the 3:1 structure deforms smoothly into a state with 
larger than half polarization by small rotations of the 
spins. Because this deformation is smooth, we expect 
that the space group symmetry of the "canted ferrimag- 
netic" state just above the plateau will be at least as 
low as the PAs32 symmetry of the R-state. This should 
be observable in neutron scattering as the persistence of 
magnetic scattering peaks - present in P4a32 but not the 
Fd3m space group of the ideal spinel - in the field region 
above the plateau. This is a distinct prediction of the 
Einstein model. 

In addition to this persistent discrete symmetry break- 
ing, this region also exhibits XY long range order of the 
spin components perpendicular to the field. The nature 
of this XY order is not apparent from simple arguments. 
Classically, it can be analyzed for the Einstein model by 
assuming small deformations of the spins, 

Si = (BeikMuny/l - IV'il 2 ) , (17) 

where crj = ±1 is the Ising magnetization of spin i in 
the R-state, and tpi is the transverse spin represented 
as a complex vector. By inserting this into Eq. i|13|) and 
expanding to quadratic order, one obtains a quadratic 
form in ipi. This can be diagonalized using Bloch's the- 
orem to obtain a set of 16 (one per site of the R-state 
unit cell) "bands" . The first vanishing eigenvalue(s) of 
this spectrum on increasing field signals the upper edge of 
the plateau. The eigenfunction (wavevector(s) and wave- 
functions) is the order parameter of the XY magnetism 
in the canted ferrimagnet. 

We have carried out this calculation for b = 0.1 and 
b' = 0.0,0.01,0.025,0.05 (the same values for which we 
plotted the magnetization curve in Fig. , and a range 
of magnetic fields sweeping through the transition off the 
plateau. Our findings are that the excitation minimum 
is at k = (the T point), with an eigenfunction which 
retains all the point group symmetries of the R-state. 
The resulting non- vanishing XY components of the spins 
just above the plateau are equal on all minority sites and 
equal on all majority sites. However, their direction is 
opposite on the two types of sites, and the magnitude of 
the transverse component on the minority sites is always 
larger than that of the majority sites, by a factor which 
is very close but not precisely equal to 3. This indicates 
that there is a non-vanishing total XY magnetization - 
transverse to the field axis. It is very small, but non-zero. 
The magnetic state is demonstrated in Fig- EJ Physically, 
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FIG. 6: (Color online) Spin configuration on one tetrahedron 
for magnetic fields just above the half magnetization plateau. 



this would be visible in a single-domain sample as a spon- 
taneously deviation of the magnetization axis from that 
of the applied field upon leaving the plateau. 

The reduced space group symmetry (relative to Fd3m) 
of this "canted ferrimagnetic" phase implies that there 
must be at least one phase transition between this state 
and the fully saturated ferromagnetic state, even at 
T > 0. Moreover, if there is only a single phase tran- 
sition, a standard Landau analysis - see Appendix [S] - 
predicts it cannot be continuous. This is not the case 
in the BP model, where for small b one observes a single 
continuous transition into the ferromagnetic state. There 
are two possibilities: a first order transition from the 
canted ferrimagnet to the fully polarized state, or an in- 
termediate phase transition between the P4^32 (R-state 
structure) canted ferrimagnet and a partially-polarized 
canted ferromagnet with higher space group symmetry. 
The latter transition is a possible explanation of the ob- 
served magnetization feature at H « 37T in HgC^O^ 4 
Some theoretical aspects of the magnetic behavior near 
saturation fields have been investigated in Ref.[23|. 



IV. FURTHER NEIGHBOR INTERACTIONS 

In this section we consider how the degeneracy on the 
plateau might be broken by further neighbor interactions. 
We take as our model the BP Hamiltonian (Eq.©) plus 
additional second and third neighbor exchange interac- 
tions: 



n 



njBP 
H-Ctt 



■h 



E 



• S-i 



j 3 Sk 

({(he))) 



(18) 



The corresponding pairs of sites were indicated in Fig|2 
Since the ground states of the BP model are exactly 
the set of 3:1 collinear states, we may treat the addi- 
tional small J2 and J3 exchange couplings in Ea. (|18l) as 



perturbations. To leading order, this amounts to simply 
replacing Sj = cr^z with Ising spins cr* = ±1 satisfying 
the 3:1 constraint. Doing so, one obtains the same effec- 
tive Ising Hamiltonian as we found in the spin-lattice 
model, Eci. (|14|) . but with the effective Ising couplings 
replaced by the physical second and third neighbor ex- 
changes, J2/3 = ^2/3! Thus it is immediately obvious 
that if J2 and J3 are chosen in a way to match those of 
the Einstein model (i.e. J3 = 2 J2), the R-state will again 
be favored. 

What state is selected by more general exchange in- 
teractions? Consider two adjacent tetrahedra. If the mi- 
nority spin is at the site joining the two, this is the only 
configuration the constraint will allow (Fig. 2(a) I . Oth- 
erwise, two minority sites must be present, in any config- 
uration. This pair can either be next-nearest neighbors 
(Fig. [2(c)) or 3 rd -nearest neighbors (Fig. |2(b]]). The lat- 



ter is preferred when J3 > J2, and otherwise the former 
tends to be preferred. This can easily be seen by calcu- 
lating the total Ising energy of both configurations. The 
Ising Hamiltonian is satisfied by having all pairs of tetra- 
hedra in the preferred ( "bent" ) generic configuration of 
the two. This is exactly the same "bending" rule found in 
the previous section, arising from coupling to lattice dis- 
tortions. Therefore, the R-state will be the ground state 
if and only if J3 > J2- If the opposite relation holds, 
J 3 < J21 "unbent" configurations are instead favored. 
This leads to a rather different state with Rim symme- 
try. In this state, the magnetic unit cell is not enlarged 
relative to the non-magnetic one. Instead, all equivalent 
(i.e. "up" or "down") tetrahedra have the same specific 
3:1 configuration. This state is thus directly analogous to 
the "uud" state expected for a kagome antiferromagnet 
in a field. It is only 4-fold degenerate (compared to the 
8-fold degeneracy of the R-state) . 

One may also study the further-neighbor exchange 
model at higher and lower magnetic fields. By straight- 
forward application of the same methods used in 
Sec lIII Bl we find that for J3 > J2, the further- neighbor 
exchange interactions favor exactly the same canted fer- 
rimagnetic state as the Einstein model just above the 
magnetization plateau (see Fig.^. 

At zero field, the situation is more interesting. Con- 
sider a pair of adjacent tetrahedra, as in Fig. [3J This 
is the smallest cluster of tetrahedra that includes spins 
interacting via further neighbor exchange. If the cou- 
pling b is assumed small relative to J2 , J3 , and J2 > J3 
then the ground state is collinear, preferring a configu- 
ration of parallel spins on a straight line, as opposed to 
a "bending ice" configuration. The J2 < J3 case is how- 
ever more complicated, and the outcome unclear from 
our limited analysis of a single pair of tetrahedra. Setting 
6 = 0, and assuming the magnetization on each tetrahe- 
dron vanishes, if J2 — J3 we find the further neighbor in- 
teractions sum to a constant. As a result we can subtract 
J2 from J3 and deal only with the third neighbor interac- 
tion terms, with the coefficient J3 — J2 > 0. Minimizing 
the third neighbor interaction terms on this single pair 
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of tetrahedra we find a non-collinear state with a mini- 
mum energy of — ^ (J3 — J 2)- The "bending ice" collincar 
state gives an energy of — (J3 — J2), significantly higher. 
In the non-collinear state, the four spins on each of the 
two tetrahedra point in four different directions defining a 
tetrahedron in spin space, i.e. with all neighboring pairs 
of spins at w 109° angles to one another. However, we 
cannot tile the entire lattice with this configuration on 
each pair of tetrahedra. Given the large energy difference 
between the collinear bending state and the non-collinear 
state on the single pair, it is reasonable to expect that 
some other non-collinear state realizes the energy mini- 
mum. Of course, as b is increased, collinear states must 
be favored even at zero field. For sufficiently large b, the 
ground states for J3 > J2 are the identical "bending ice" 
configurations found in the Einstein model. 

Other possible interactions in the chromites include 
dipolar interactions. Spin ice compounds, for exam- 
ple, often have appreciable dipole-dipole interactions^ 
that cannot be ignored. Using the estimate of the mag- 
netic moment of Cr +3 ions in CdCr 2 04, fi ~ 3.7^s 
we estimate the dipole-dipole coupling strength to be 
D = j^^r- - where r nn = 8.567A is the near- 

est neighbor distance in the CdCr204 pyrochlore struc- 
ture. With the experimental observation of the magnetic 
ordering temperature at a few degrees Kelvin, it is ev- 
ident that the energy scale of the dominant degeneracy 
breaking mechanism, regardless of its details, is at least 
an order of magnitude larger than that of dipolar in- 
teractions. For this reason we believe that the dipolar 
interactions play no significant role in the physics of the 
chromite spinels, and can be safely neglected. 

V. DISCUSSION 

In this paper we have studied classical mechanisms 
of degeneracy breaking in pyrochlore anti-ferromagnets 
and how this may be related to certain features of 
the low temperature magnetization curves of the spinel 
chromites, ZnCr204, CdCr204, and HgCr204. A 
very simple Einstein phonon model predicts a half- 
magnetization plateau with a unique ground state spin 
configuration. Given the known presence of large spin- 
lattice coupling in these materials;^ this seems the most 
likely explanation of the plateau state. Neutron scatter- 
ing measurement of the R-state structure on the plateau 
in any of these materials would provide good support of 
this proposal. 

Nevertheless, the same ground state could in princi- 
ple be stabilized by other sorts of interactions. For in- 
stance, as we have shown, further neighbor exchange with 
J3 > J2 would also lead to the R-state. This relation 
amongst exchanges seems contrary to simple expecta- 
tions that antiferromagnetic superexchange decays with 
distance. However, the exchange paths for these two in- 
teractions (and indeed some further-neighbor exchanges) 
are quite similar (see Ref ill , so one should have an open 



mind to this possibility. It would thus be desirable to 
have an independent comparison of the two theoretical 
models. 

One further check on the Einstein model would be to 
consider its predictions in zero field. As we have seen, 
despite the selection of a unique ground state on the 
plateau, in zero field the Einstein model continues to pre- 
dict a large degeneracy of states. In reality, of course, 
this degeneracy will be broken by further interactions 
(e.g more complex phonons, further-neighbor exchange, 
Dzyaloshinskii-Moriya interactions) beyond the Einstein 
model. Interestingly, the zero field ordered phases for 
ZnCr204, CdCr2 04, and HgCi'204 are known from neu- 
tron scattering and are known to be different&Z^ This 
is indeed consistent with the Einstein model in the sense 
that the further very small interactions beyond the model 
would be expected to select different states in each ma- 
terial. Interestingly, the Einstein model suggests that 
despite this panoply of phases at zero field, all these ma- 
terials may display a universal ordering on their magne- 
tization plateau: the R-state structured 

As we have seen in Sec lIVI under some circumstances 
(J3 > J2 and b not too large), the further- neighbor ex- 
change model may have a non-collinear ground state. 
This would be clearly distinguishable from the collincar 
states preferred by the spin-lattice interactions. For 
instance, Goldstone's theorem implies that a collinear 
ground state in zero field will have two gapless spin 
wave modes, while a non-collinear state will have three. 
However, we caution that additional effects beyond the 
Einstein model, especially Dzyaloshinskii-Moriya inter- 
actions, could induce some small non-collinearity even if 
the predominant interactions are of the spin-lattice type. 

Another general prediction of the Einstein phonon 
model is that the interactions (b) which stabilize the 
plateau are of the same order as those (b') which select 
the ordered R-state out of all possible plateau configura- 
tions. This has physical consequences. Specifically, the 
energy cost of a spin excitation which leads to a devi- 
ation from the quantized plateau magnetization is ex- 
pected to be of the same order as a spinless excitation 
which re-arranges the 3:1 configurations but leaves the 
magnetization unchanged. On heating the sample, the 
former excitations are responsible for the rounding of the 
plateau, while the latter are responsible for the destruc- 
tion of the R-state magnetic order, i.e the restoration of 
the spinel space group symmetry. Because both excita- 
tions will be excited roughly equally, we expect that the 
thermal phase transition from the R-state to the high- 
temperature phase should occur at a critical temperature 
T c which is of the same order of magnitude as the scale 
at which the plateau forms, and the 3:1 constraint itself 
is rather strongly violated. 

Because the 3:1 constraint is not significant at this tem- 
perature, a conventional Landau-Ginzburg- Wilson anal- 
ysis of this critical point is valid (see Ref.[0| for a dis- 
cussion of the alternative scenario which would apply 
if T c <C T p , where T p is the temperature at which the 



11 



plateau forms) . The result of such an analysis - detailed 
in Appendix ^ - is that the thermal transition should 
be first order. This is in agreement with experimental 
findings. 

There are many open directions for future work. An 
important one is to understand more microscopically 
the mechanisms of exchange interactions and spin-lattice 
coupling. Goodenough-Kanamori analysis actually pre- 
dicts a competition between two processes effecting the 
nearest neighbor exchange interaction: antifcrromagnetic 
direct Cr-Cr exchange, and ferromagnetic Cr-X-Cr su- 
perexchange in the ideal crystal structure, because of 
90° Cr-X-Cr (X=0,S,Se) bonds^ The angle of the Cr- 
X-Cr bonds is generally not exactly 90° however, and 
is affected by the it-parameter in the spinel structure. 
The prevailing belief is that the competition of these 
two processes is most strongly effected by the overall ex- 
pansion/contraction of the lattice, but the different u- 
parameters in different materials may also be important. 
Longer distance super-exchange processes between sec- 
ond, third, and fourth neighbor pyrochlore sites involve 
comparable exchange paths, and their relative magni- 
tudes are not presently clear. Turning to spin-lattice 
coupling, an interesting speculation is that the impor- 
tant phonon modes are those which modify the Cr-X-Cr 
angles, thus strongly affecting the nearest-neighbor su- 
perexchange contribution. This is an appealing possibil- 
ity given the observation of changes between ferromag- 
netic and antiferromagnetic behavior in small changes 
of temperature and field in HgCr2S4. ? - A broader un- 
derstanding of the microscopic mechanisms of spin and 
spin-lattice interactions in these materials would be of 
considerable interest even beyond the realm of frustrated 
magnetism. 

In summary, the study here clearly highlights the sen- 
sitivity of the magnetic state in antiferromagnetic py- 
rochlores to further microscopic interactions. A simple 
Einstein model makes the prediction that the R-state or- 
der should be observed in many materials and we hope 
that detailed neutron scattering studies will be forthcom- 
ing to test this prediction. 
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APPENDIX A: LANDAU-GINZBURG THEORY 

FOR THE FINITE TEMPERATURE 
TRANSITION OUT OF THE PLATEAU STATE 



In the R-state^ the enlarged spin periodicity, includ- 
ing 4 unit cells of the underlying pyrochlore lattice, man- 
ifests itself in the appearance of non-zero Fourier compo- 
nents with a momentum vector not at the Y point in the 
Brillouin Zone (BZ). These Fourier components there- 
fore serve as order parameters for the R-state. Then, 
using these order parameters and their transformation 
under the various lattice symmetries, we can construct 
a Landau-Ginzburg theory to determine whether a finite 
temperature transition out of this ordered state should 
be 1st order or 2nd order, at least at the mean field level. 

Previous work by some of the authors^ has encoun- 
tered the same R-state in a different formulation, using 
a PSG (Projective Symmetry Group) analysis. We find 
the Fourier components are at the 3 momentum vectors 
k a = §(7r,0,0), k 2 = i(0,7r,0), and k 3 = ±(0,0,7r). 
These form a k-star of the point symmetry group of the 
lattice. This particular k-star has only 6-dimensional ir- 
reps (irreducible representations) of the point symmetry 
group. Therefore the order parameter can be cast as a 
6-component real vector. These components can be un- 
derstood as degrees of freedom equivalent to the real and 
imaginary parts of the 3 Fourier components. 

The R-state corresponds to 8 configurations of this 
more general order parameter v = (<ti, 02, 03, 0, 0, 0) v 
or v — (0, 0, 0, tJ4, 05, ae) v where <jj = ±1, and we allow 
only the 8 cases where Cio^o^ = +1 or o^a^a^ = +1. 

This order parameter was derived using the PSG anal- 
ysis. However, it is possible to derive this order param- 
eter more directly, by finding the "density" of minority 
sites pj = I (l — Sj). A straightforward Fourier analysis 
results in 



p(x,y,z) = 

— (l + (T4V2 COS — X + (T5V2 COS — IJ + <7g V2 cos — z 

4 V a a a 

— (ii\/2sin —x — er 2 -\/2sin — y — a^V^sin —z) , 

a ' a a ) 

(Al) 
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where (a;, y, z) are the real space coordinates of the py- 
rochlore lattice sites. 

In this representation we construct invariants (or 
Casimir operators) from the 6 components {vj}® =1 of v. 
The Landau-Ginzburg (LG) theory is then constructed 
out of all these invariants. In this way we find the most 
general LG free energy allowed by the symmetry of the 
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lattice 



F —mv + 7 (V1V2V3 + V4V5VQ) + u\v 



-4 



3 6 
j=l t=4 



(A2) 
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The 8 R-state configurations are favored for a part of the 
coupling space 7 < 0, m < 0, u\ > 0, u 2 > 0, and u 3 < 0. 

The transition into the orderless state is tuned by m 
changing sign. We take v — (1,1,1, 0, 0, 0) v without loss 
of generality, and simplify the free energy to a single vari- 
able function 



'L.4 



F(v) = mV - |7|iT + \u'\v 



(A3) 



Due to the cubic term, the transition is predicted by MFT 
to be 1st order. 



APPENDIX B: PROOF OF BENDING RULE 
INDUCING THE R-STATE 

In this appendix we explain how from the "bending" 
rule for the 3 : 1 configurations, we can construct only 
the R state. 

Consider the unit cell of the R state, including 4 up 
pointing tetrahedra as in Fig-QJ We will show that using 
bending configurations on pairs of tetrahedra, we will 
find a unique state (up to lattice symmetries). 

We start by picking a pair of tetrahedra to be in the 
bending configuration. We can then pick the unit cell 
orientation such that it m atches the placements of the 

the two adjacent tetrahedra 



minority sites in Fig. 7(a) 



share site 1 and the two minority sites are 2 and 6. 

Next we consider the tetrahedron pairs sharing sites 3 
and 5, shown in Fig. 7(b) For the pair sharing site 3, we 
can pick sites 4 or 7 to be minority sites. Similarly, for the 
pair sharing site 5, we can pick sites 4 or 8 to be minority 
sites. We cannot pick 4 and 7 or 4 and 8, since then we 
will have nearest neighbor minority sites. We also cannot 
pick 7 and 8, since then the pair of tetrahedra sharing site 
4 will not be in a bending configuration - sites 7,4 and 8 
sit on a straight line. We must therefore choose site 4 to 
be a minority site, (see Fig. |7(c)] |. Already a tendency to 
form flippable plaquettes is evident. 

Now we turn to the upper layer of tetrahedra (the 
tetrahedra outlined by dashed lines). Considering the 
pair of tetrahedra adjoined at site 7, we conclude that site 



9 cannot be a minority site. Similar considerations deem 
sites 10 and 11 cannot be minority sites. This should be 
most clearly evident from the 3-fold rotational symmetry 
of the minority site configuration about an axis perpen- 
dicular to the paper. We mark these sites in Fig. |7(d)| by 
a full gray circle. 

The pair of tetrahedra sharing site 7 can allow a mi- 
nority site o n either site 12 or 13. If we choose site 13 (as 
in Fig. 7(e) I, then now we cannot choose site 14 to be a 
minority site, as it neighbors site 13, and we also cannot 
choose site 15 to be a minority site, as 13,14 and 15 sit 
on a straight line. Now considering the pair of tetrahedra 
sharing site 16, we reach an impasse - the tetrahedron of 
sites 11, 14, 15 and 16 cannot have a minority site on any 
one of its corners! Therefore, we cannot choose choose 
13 to be a minority site, we can only choose site 12! 

Due to the 3-fold rotation symmetry of the spin struc- 
ture we have layed out so far the same argument ap- 
plies to all three down pointing tetrahedra in the upper 
(dashed) layer. The single up pointing tetrahedron must 
therefore have a minority site locate d at s ite 17. The re- 
sulting configuration, shown in Fig. 7(f) is the unit cell 
of the R-state. 

Next we demonstrate that the R state is the only 3 : 1 
spin configuration that can be constructed from these 
cell structures. We can add the minority site positions 
on the remaining tetrahedra in the cell cluster, as shown 
in Fig. 18(a)] Each cell has 4 hexagonal plaquettes. One is 
"flippable" - alternating between majority and minority 
sites. The other three plaquettes have one minority site, 
and 5 majority sites. There is only one way to attach 
an identical cell with the flippable plaquette, but apriori 
there are three ways to attach two cells through a pla- 
quette with only 1 minority site (we have the freedom to 
choose which one of the other 3 plaquettes of the second 
cell is the flippable one). 

Consider the cell sharing the hexagonal plaquette that 
involves sites 1,6,16 and 14. Since sites 6 and 15 are 
already set to be minority sites, the flippable plaquette 
in this cell must be the one going through sites 6, 16 and 
15, marked by purple lines. We therefore do not have 
three choices of how to place the minority sites on this 
cell, but rather only one. The whole cluster, with the 
given spin configuration, still retains a 3-fold rotational 
symmetry about an axis perpendicular to the paper. We 
use this symmetry to deduce that the same considerations 
apply to the other three plaquettes with one minority 
site, and there is no freedom in how to choose the spin 
configurations in all the surrounding cells. Therefore, 
there can only be one way to arrange the minority sites 
with the given cell, and that is the R-state configuration. 
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FIG. 7: (Color online) Stages in the graphical proof that the 
bending rule allows only the R state to form with the 3 : 1 
constraint. The down pointing triangles, with lines joining at 
their centers, represent tetrahedra pointing out of the plane 
of the paper, with the center point representing the upper 
corner. The up pointing triangles in the figures represent 
tetrahedra pointing into the paper, with the lowest corner 
being represented by the center position (sometimes marked 
by a small cross). The dashed lines represent an upper layer 
of tetrahedra, above the solid line tetrahedra. It is instructive 
to contrast this projection with Fig. Q- both show the same 
cluster of tetrahedra. To avoid clutter, we mark only the 
minority sites, by down pointing (red) arrows. We mark sites 
that we conclude are not allowed to be minority sites by a full 
gray circle. 
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FIG. 8: (Color online) Stages in the graphical proof that the 
bending rule allows only the R state to form with the 3 : 1 
constraint. The down pointing triangles, with lines joining at 
their centers, represent tetrahedra pointing out of the plane 
of the paper, with the center point representing the upper 
corner. The up pointing triangles in the figures represent 
tetrahedra pointing into the paper, with the lowest corner 
being represented by the center position (sometimes marked 
by a small cross). The dashed lines represent an upper layer 
of tetrahedra, above the solid line tetrahedra. It is instructive 
to contrast this projection with Fig. 0- both show the same 
cluster of tetrahedra. To avoid clutter, we mark only the 
minority sites, by down pointing (red) arrows. We mark sites 
that we conclude are not allowed to be minority sites by a full 
gray circle. 



